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A theory is presented for the modification of the line-shape functions and absorption coefficient due 
to the breakdown of the quasistatic approximation. This breakdown arises from the effects of 
molecular motion and increases the absorption in the near wings. Numerical calculations for the 
high-frequency wing of the v 3 band of C0 2 broadened by Ar are reported and it is shown that these 
effects are significant near the bandhead. The importance of such corrections in other spectral 
regions and for other systems is discussed briefly. 




I. INTRODUCTION 

The calculation of accurate frequency profiles for mo- 
lecular vibration-rotational transitions valid over a range of 
different physical parameters (e.g., number density, tempera- 
ture, etc.) remains an unsolved problem. Various theories 
valid for small frequency displacements from the resonant 
line center have been proposed over the years. 1 For instance, 
the impact theory in which the duration of collisions is as- 
sumed to be zero leads to Lorentzian-shaped profiles that 
apply quite generally to non-overlapping lines arising in low 
density gases. 2 ' 3 On the other hand, only relatively recently 
have tractable theories been proposed that are valid for large 
frequency displacements. Rosenkranz, 4 ^ starting the basic 
formalism given by Fano, 6 has developed a far-wing theory 
based on the quasistatic and binary collision approximations 
and was able to obtain theoretical results for self- and 
N 2 -broadened H 2 0 spectra (the so-called continuum absorp- 
tion) which are in substantial agreement with the experimen- 
tal data. 7 But because of the assumptions made, this theory is 
applicable only to the high-frequency wing of the pure rota- 
tional band. Rosenkranz’s theory has been subsequently ex- 
tended and improved by other authors. 8-13 Boulet et al. s pro- 
posed a line-by-line generalization of Rosenkranz’s theory, 
but because of the extensive computations involved this ap- 
proach has been implemented only for atom-molecule 
systems. 14 Ma and Tipping 9-13 in a series of papers have 
considered a number of improvements, removed the restric- 
tion for the applicability only to high-frequency wings, and 
applied their results to the calculation of self- and foreign- 
broadened H 2 0 for frequencies up to 10 000 cm -1 . Their 
results are generally in good agreement with both laboratory 
and atmospheric data. 9-13,15 However, in all the formalisms 
mentioned above, the time dependence of the interaction be- 
tween an absorber and a perturber is ignored. • 

In the present paper we consider corrections to the qua- 
sistatic theory implicit in all the theoretical treatments dis- 
cussed above. 4,5,8-15 In Sec. II, we present the general for- 
malism and give results for the effects of molecular motion 
(neglected in the quasistatic approximation) on the spectral 
density and on the line shapes. This formalism is applied in 


Sec. Ill to the case of the band of C0 2 broadened by Ar. A 
discussion of the numerical results and conclusions of the 
present work are presented on Sec. IV. 

II. THEORY 

A. The basic approximations 

In order to take into account the molecular motion for 
the wing region, we make some basic approximations be- 
sides the binary collision approximation. 

(1) The quasistatic approximation. In the quasistatic ap- 
proximation, the duration of a collision is assumed to be 
infinite and the positions of interacting molecules are fixed in 
space during the period of time of interest. One usually treats 
the relative motion of the two molecules classically; i.e., the 
ensemble average over the relative motion is carried out by 
integrating over r, the separation between two interacting 
molecules, with a statistical weight for a pair of molecules 
p('m) (j.} related to the interaction between them. We note that 
since the interaction is independent of time, the molecular 
motion is totally ignored. This approximation limits the ap- 
plicability of the theory to the far-wing regions. In the near- 
wing region, although we can still assume the duration of the 
collision is very large, we cannot take it as infinite as we did 
previously. Therefore, the molecular motion should be taken 
into account during the time of interest. 

We assume the duration of collision, r, is very large. 
Therefore, during the time of interest (f~l/Ao>, which is 
small), the displacements of molecules are small and the in- 
teraction between them changes only slightly. We are thus 
able to write the interaction as two components: a main com- 
ponent V independent of time and a small correction term 
V,(t) that is dependent on time. Consequently, the total 
Hamiltonian of two interacting molecules similarly consists 
of two parts, a main part independent of time and a small 
time-dependent correction term 

H(t)=H a +H b +V+V l (t)=H+V l (t); (1) 

here the subscripts a and b refer to the absorber and per- 
turber, respectively. 
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(2) The closest region approximation. The interaction be- 
tween two molecules depends on their separation. In general, 
during the whole collision process, the closer the distance 
between two interacting molecules becomes, the stronger the 
anisotropic interaction. As is well known, the anisotropic in- 
teraction is mainly responsible for collision broadening 
which is related to the coupling between the interaction and 
the transition of the absorber molecule during the collision 
process. Therefore, the strongest coupling effect takes place 
when two interacting molecules are near the separation of 
closest approach (we choose that time as t=0). In compari- 
son with other parts of the collision trajectory, the contribu- 
tion of this coupling comes mainly from the closest region, 
and we can assume that it can be ignored elsewhere. 

(3) The classical ensemble average approximation. We 
do not treat the relative translational motion of two interact- 
ing molecules quantum mechanically, rather we consider 
their motion as a classical trajectory. In this case, the en- 
semble average over the relative motion is an ensemble av- 
erage over all classical trajectories. However, based on the 
closest region approximation, the coupling effects from the 
whole trajectory arise mainly from separations around the 
closest point r 0 . We introduce a parameter ij to estimate the 
fraction of the trajectory consistent with the closest region 
approximation that contributes significantly. With a straight 
line trajectory approximation, one is able to carry out this 
ensemble average. Based on the approximations discussed 
above, we can develop a formalism to take into account the 
effects of the molecular motion on the line wing shapes. 

B. General formalism 

The correlation function of the dipole moment is defined 
by 


From Eq. (1) and the perturbational method, we obtain 

t/(r) = e' 7, ' / *Jl + t|V I (r , )dt' 

Jo Jo 


where 


( 6 ). 


(7) 


In the above expression, we have introduced the Liouville 
operator L defined by 


L = {HI* —IH*)/h. 
Then, we have 

= 1+i 


(8) 


MO) 


i\‘v l l (t')dt' + 

Jo 

if , V[(t')dt' + ■ 
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= e ,x X 0 ) + ie iH, \j V, ( t ' ) dt ' , 0 ) 


x 1-t 


e~ iHl + ■ 


=/*o (t) + ie iLl \'v[(t')dt',rtO) 
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where 

ft o (t)^e iH 'tt(0)e- iH '=e iL 'M(0). 
It is obvious that 


(9) 


(10) 


"V)>=2 p,<<U" , (o)-/4 f ' ) (/)|;>, (2) 

i 

where p, is the normalized Boltzmann factor of the initial 
state i. In the above expression, is the dipole opera- 

tor of the absorber molecule in the Heisenberg representation 
which satisfies the equation 

= (3) 

where the square brackets indicate commutating operator. In 
the present case, the Hamiltonian H(t) explicitly depends on 
time and we are unable to simply express f4, H \ t) as 
e' w (')' M p4 w )(0)e~' w (' ,,M as we did previously. However, we 
are able to solve Eq. (3) using a perturbational method as 
follows. For simplicity, we will omit the subscript a of 
ft indicating the absorber molecule and the superscript 
(H) indicating the Heisenberg representation. We introduce a 
time displacement operator U{t) and assume that 

M (t) = C/(t)M(0)t/(t) t . (4) 

With Eqs. (3) and (4), we obtain 

U{t) = iH{t)U{t)/h. (5) 


/«o(0) = /40). (11) 

The correlation function can be written as 

Tr{pA*(0) ■ A*(0} = Tr{pA^)(0) • e ,u ^o(0)} + Tr| p/u 0 (0) 

.ie iL \[‘v\{t')dt'\^{0) 

[Jo 

+ -•, (12) 

and its Fourier transform, the spectral density F( o>), as 

F(oj) = -— P e io "Tr{pMO) M(0}^ = F 0 (a.) 

Zir j - ® 

+ Fi(w)-l . (13) 

The zero-order term of the spectral density, F 0 {w), is defined 
by 

Fo(w) = 2^ S-J ia " Tr ^ (0) - e ' X 'Fo(0 )}dt 

= ^ Re L e ‘ a,Tl ^ e ~ iU PM}dt, (14) 
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where for simplifying the notation, we have omitted the sub- 
script and the zero-time argument of the dipole moment op- 
erator /io(0). The first-order term of the spectral density can 
be expressed as (see Appendix A) 







dt. 



(15) 

where V'j(r) is a Liouville operator corresponding to Vj (/). 

It is well known 9-13 that since the isotropic interaction 
does not depend on the rotational quantum numbers, the 
Liouville operator associated with it will cancel out when it 
acts on the line space. Therefore, the isotropic interaction 
need not be included in the Liouville operator L nor in its 
time-dependent correction part Vj(r). 


C. The interaction at the closest region 

In the present paper, we assume that the anisotropic in- 
teraction takes a product form 

V'aniso[r(f)]=^[r(t)]G an i SO [0(t), ¥>(/)], (16) 

where r(t), 6(t), and <p(t) are the spherical coordinate com- 
ponents of r(r). We note that for simplifying the notation, the 
orientations of the two interacting molecules have not been 
explicitly indicated in the above expression. Since the vector 
r (r) is time dependent, the anisotropic interaction explicitly 
depends on time. However, based on the closest region and 
the quasistatic approximations, we are only interest in the 
behavior of this interaction around the closest point r 0 and 
for time f<Sr. We choose r 0 as the k axis of a space-fixed 
frame (i, j, and k) and r (t) is on the plane with <p=0. We can 
expand V anjso (r,d) around f=0 as 


j dr d d$ d\ , 1 

o (r,0) = V' ani Jr o ,O)4-r|— — +— — j V' miso (r,0)| /=O +- t 

v ™*o{r,0)\ t=0 +--- 


dr\ 2 d 2 l dr )l d0 \ £ d 2 r d 

dt I dr* + dt ) 1 dt) drdd + dt 2 dr 


ldd\ 

2 d 2 

d 2 e 

d 

1 dt) 

dd 1 

+ d? 

Je 


dd d 

~ l^aniso( r 0 1 0 ) l^aniso( r 0 > 0 ) 


4 ' ! 


d 2 r d ldd\ 2 d 2 d 2 0 d 


dt 2 dr + \d7) Je* + dFde 




(17) 


r 


where in the last step of deriving Eq. (17) we have used the 
fact that ( dr/dt)=0 at the closest point. Since the set of 
trajectories has rotational symmetry with respect to the axis 
k, we expect that the contribution from the second term will 
be zero and we will simply omit this term in Eq. (17). Con- 
sequently, we are able to express V miso (r,S) as two parts: a 
dominant part ' / a niso( r o.°) which is independent of the time 
and a small correction part which is proportional to t 2 . It is 
worth mentioning that the correction term is an even function 
of time. Thus, we rewrite Kj(f) as 

V\(t) = t 2 V \ , (18) 


where 


dr dd d 2 d 
m\ 2 — — +r “tt 
dt dt dr 


1 d 

TJe 


_ ja a l^ani so(^>^)> 


( 21 ) 


where m is the reduced mass of the two interacting mol- 
ecules. Then, we obtain 


d 2 r 

d?\ 


r = r Q , 9=0 


I dd\ 2 Id Id 
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1,2 2 
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dt) dd 2 
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d 2 d d 

Je 
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(19) 


Based on classical dynamics, we have 

d 2 r ldd\ 2 d d 

m 1? =mr [ ~dt) ~d? V *° < ' r) ~d? V ^ r,9 ‘ ^ 20) 


In the last step of deriving Eq. (22), conservation of the 
energy 

1 7 (dd\ 2 1 , 

— mrgi dt I — 2 m V ® — ^iso( r 0 ) ^aniso( r 0 >0), (23) 

has been taken into account. In the above expression, i/ x is 
an initial relative velocity when the separation between two 
interacting molecules is infinity. Since \mrl(dd/dt) 2 3=0, the 
range of r 0 which can be approached is limited. This fact will 
be used later. From Eq. (21) we obtain 


J. Chem. Phys., Vol. 100, No. 8, 15 April 1994 


5570 


Q. Ma and R. H. Tipping: Quasistatic line-shape theory 


<f_d\ 

Tt 


r=r a ,e= 0 


fltTg 36 


^aniso(^0>®)- 


Using Eqs. (22) and (24), V\ can be expressed as 


Vi- 


2 mr n 


mv x -\2 + r 0 — Jt^rg) 


(24) 


"*■ ^aniso( r O>0)]J ^aniso ( r 0 > 0 ) 

1 l dd V , 1 , 

+ 2 (dr) ° Vini ” 2mr\ - (25) 

where 0 2 V^3 2 V/dd 2 and (OV) 2 ^(3V/36) 2 . We note that 
the first term of Eq. (25) depends on the velocity v x while 
the last term is related to a torque. 16 For simplifying the 
notation, we will omit the subscript of r 0 in the following 
discussion. 


In the above expressions we explicitly denote the summation 
indices of the trace, a or b, for the absorber or bath molecu- 
lar variables, respectively. 

Now we focus on the correction term F^wi). With Eqs. 
(A4), (AS) and (18), we have 


1 f« 

~i(co) = — Re Tr 
tt Jo 

»- Re Tr f 
tr Jo 


H ■t 2 e i{ “~ L), V[ — pfx 

(i) ~~~ L 


ft ■rV ( "" i o)'V r ; - — 7 PH 


,-L 


dt, 

dt, 

(28) 


where the time-independent Liouville operator L consists of 
two components, 


L =L 0 + ^i • 


(29) 


D. The effect of the molecular motion on the 
expression of F,( «) 

In a previous paper, 10 we have derived an expression for 
F 0 ((o) in which, at least formally, the bath variables have 
been isolated 


1 

F 0 (co) = — Im 



dco' Tn 


F 


1 

co— co' —L a 


corresponding to the unperturbed Hamiltonian of the ab- 
sorber and bath molecule and the time-independent interac- 
tion, V^ro.O), respectively. We assume that 

(30) 

where the limitation of this approximation has been dis- 
cussed in a previous paper. 10 Based on it, we can rewrite the 
resolvent operator [!/(«—£.)] in Eq. (28) as 
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where 
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(32) 


In order to separate the summations over the absorber vari- 
ables and the bath variables, we first explicitly write the sum- 
mation indices of the trace in Eq. (32) 


— Tr 


y' 

( co-LqY 1 co'-Li 'co-co'-Z,, 
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PM 


Then we can carry out the summation over the bath variables 
of the trace using the fact that there is no contribution from 


J. Chem. Phys., Vol. 100, No. 8, 15 April 1994 


Q. Ma and R. H. Tipping: Quasistatic line-shape theory 


5571 


the Liouville operator whenever it acts directly on ft. 
Therefore, the trace in Eq. (32) can be written as 


Tr 


where 


(oj- 


~L J 7 V 1 u'-lJ. o-o'- 


P<‘V 


(33) 


o'-L x 


v[ -jzr) ~Tr| V[ —r~r P 




o’ -Li 


(34) 


Thus we obtain an expression of Fj(o>) in which the bath 
variables have also been isolated just as in the expression of 
F 0 M 


Gl «)03| = (G a - G $ )\ a)(/3\ a><0| . (38) 

Then, the matrix elements of the Liouville resolvent operator 
in the representation constructed from its eigenvectors can be 
written as 


^dzr 1 (39) 


Based on the knowledge about the resolvent operator 
[l/(c»> — Z. !)] obtained above, we now discuss the resolvent 
operator 
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2 1 

F,(o»)« Im — do' Tr 

it 2iri J -oo 
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X w\ 


1 o'—Lij o—o'—L, 


P (a V 


(35) 


In comparing Eq. (35) with Eq. (26), we find that the mo- 
lecular motion introduces complexities in the trace over per- 
turber variables in the calculation of the spectral density. 


E. The resolvent operator (V\[M(o' - £.,)])& and its 
trajectory average 

The resolvent operator [1 /(w-LJ] is defined only in the 
representation constructed from its eigenvectors. Therefore, 
we have to find the basis in which L , is diagonal. The Liou- 
ville operator L ( corresponding to the time-independent part 


in Eq. (35). It is an operator acting on the line space of the 
absorber molecule a. Therefore, in general, it can be ex- 
pressed as 



(40) 


By explicitly writing the summation over the bath molecular 
variables, its components can be written as 




= 2 (dij\)(hil)\V'i 
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We denote by (or) and G a , the eigenvectors and eigenvalues, 
respectively, of the matrix operator G aniso ; that is 

C in Ja) = C», (37) 

or, equivalently, for the Liouville operator G in line space 
notation 


We note that the summation over the bath variables in line 
space in Eq. (41) includes only the states | i 2 il) because of 
the assumption that the bath molecule is not involved in the 
transition. 

Based on our classical ensemble average approximation, 
the ensemble average over the classical trajectories (indi- 
cated by a subscript traj.) is replaced by a weighted integral 
over r times the parameter tj. Therefore, the trajectory- 
averaged component of the resolvent operator {V r ^[l/(u> , 
- L [)])(,, traj. can be written as 
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«i, 
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VWr) 

imvl 


11/2 


is a statistical weight 17 and n b is the number density of the 
bath molecules. We note that in Eq. (42) the lower limit of 
the integration is not zero but rather the turning point r min 
which comes from conservation of energy and consequen- 
tially depends onu., the relative velocity of two interacting 


molecules at infinity. Since both the operator V[ and r min 
depend on the velocity v„, when we carry out the trajectory 
average, we have also to take into account velocity averag- 
ing, assuming a Maxwellian distribution function. However, 
for simplifying the calculations, the average over i/» will not 
be treated exactly, but approximately. We simply replace 
by a properly chosen value v„ in the corresponding expres- 
sions. Using Eqs. (25) and (39), the matrix elements 
(a/3 f |V;[l/(«/ - L x )]\a'p'') in Eq. (42) can be ex- 
pressed as follows: 


( a 0'\V[ -r^~ \ a 'f3'') = (af3'\V[\a'f3'') -r^7?r 

(x) —L J (J) -jK(r) U a 'B' 
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(43) 


where G*g = (G a + Gg ) . It turns out that the contributions to 
F x {u>) from the last two terms of Eq. (43) approximately 
cancel. In fact, if we replace |( dd/dt ) 2 by ( kT/mr 2 ) and 
then perform a trace operation, they can be approximated by 

IcT 1 
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However, in the present calculation we are interested in the 
imaginary part Im(aj8 t |V r '[l/(w"-L,)]|a'/?' t ) only. With 
the well-known formal identity 
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u~R(r)G a g+ie 
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and ((O V) 2 ) is the mean squared torque defined by 
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where P indicates the Cauchy principal value of the integral 
and <5 is a Dirac delta function, it can be expressed as 
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where ( ) n indicates the average over all orientation of two 
interacting molecules. As is well known classically 


kT(0 2 V) = ({OV) 2 ), 


(46) 


and, therefore, the contributions from these two terms cancel 
out. Based on the argument given above, we need only to 
focus on the first term of Eq. (43) and we obtain 


Xr — R{r)G a gS[(o' - R{r)G a ' g']S aa ' ■ (49) 


Finally, using Eqs. (42) and (49) we can write (see Appendix 
B for the details) 
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where 
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F. The modification of the line shapes 

Using Eq. (B7), the correction term of the spectral den- 
sity, F^oi), can be rewritten as 
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There are a lot of individual ,-,y, functions involved in 
Eq. (52). However, in practical calculations we can introduce 
intensity-weighted averaged functions to replace them. The 
intensity weighted averaged function fd(<a) is defined by 
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In analogy to our previous work on band-average line shape 
functions, 10 " 12 we can define the symmetric function 


The explicit form of the ?/((o) function is given by 
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af3 m 


(55) 


' , V-J.r) 

1/2 

dR(r ) 


1 

kmvl 

Sign 

dr ° ap 

r 


1 -- 


mv. 


d\ 


X|2 + r— j ( V iso (r) + /?(r)G^) 

x(a|p ( V“Vml/3) / 2 P/ IF,; 1 2 


'=' a pM 


(56) 


In terms of this function, the spectral density can be rewritten 
as 


F,(o*)=- 2 e~^-^ l2kT 
7r -o 

1 


1 

7 77 #*( “ “ «.;) 

(W-W.y)- 5 1 


(W+W,;) 


7 M(u>+u>ij) 


PilMijl 2 - (57) 


In the present formalism, the absorption coefficients 
a(w) consist of two parts: the main part a 0 (aj) related to 
F 0 M and a correction part o' 1 (o>) related to F^w). The main 
part of absorption coefficient obtained in the previous 
papers, 10 " 12 is given by 

4tt 2 

a 0 (<o) = - Jc ° iria tanh (^ w /2^7')[F 0 (w)+F 0 (- w)] 

^ ^ <vsh(hco/2kT) 1 
~ na J J >0 Sii ^shih^nkT) n 


1 1 . . 
lx{ „- Ui j) + 1 ——p X {u>+Ui i ) f, 


(<Ji- (Oil) 


(58) 


where 

x(<o) = 4tt 2 /i 6 w 2 2 HFml£>* 

a (3 m 

Xe -ha,[iC a +G 0 )/(,G a -G 0 )}/2kT 


2 P/IF,;| 2 . 

(53) 

V 1 

R(r) 

„-V; so (0/*r 



A i r 

dR(r)/dr 

C | 


where r a ^<d) are the positive solutions of the equation 
<o-R(r)G a p=0 which are larger than r min . In terms of 
fd(u >), Eq. (52) can approximately be expressed as 


r=r a s(") 

X(a|p (i,) p <a Vml/S) / 2 Pil F,;l 2 ’ 

/ ‘7 


(59) 
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where r a ^<d) are the positive solutions of the equation 
<t>— R(r)G a p=0. In Eq. (58), the line strength 5, ; is defined 
by the usual expression 10-12 


4tT 2 <U;; 

S '' = 3^q7 


8j\t*ij i 2 (e 


-Ej/kT_ 


e 


-EJkT} 


(60) 


On the other hand, using the results derived above the cor- 
rection term can be written in the form 


4 IT 2 

a ‘( w )~3 He 0}n “ tan h(* w/ 2Ar)[.F 1 (<i>)+F I (-a))] 


= «<. 2 




(Dsh(h(D/2kT) 1 
(i>^sh(h(o^/2kT) tt 




(<o + <*>,,) 


7 ?/j{u)+ ui,j) 


(61) 


In comparing the expressions for and a^tu), we find 
that they are similar except for the additional power of 
(oi—Wjj) and {(o+cjjj) in Eq. (61). This is consistent with the 
expectation that the correction contribution should decrease 
faster than the leading term as the frequency displacements 
increases. 


III. NUMERICAL RESULTS FOR THE Ar-BROADENED 
i/ 3 BAND OF C0 2 


As an example, we consider the effect of molecular mo- 
tion on the wing of the Ar-broadened v 3 band of C0 2 for the 
temperature T— 296 K. In this case, the bath molecule is 
simply an Ar atom and its internal quantum states need not 
be considered. In order to simplify the calculations, we 
model the isotropic potential between C0 2 and Ar by the 
Lennard-Jones model 


v 'i*>('') = 4e 




(62) 


The expression of the anisotropic potential, Eq. (16) is the 
product of the functions 


R(r)= 4e 



(63) 


and 


G ani so=^2(cos 0), (64) 

where 0 is the angle between r and the C0 2 molecular axis, 
and R 2 and A 2 are adjustable parameters. We first adopt the 
same values used by Boissoles era/. 14 with cr=3.91 A, elk 
= 153 K, /? 2 = 1.09 and A 2 = 0.266. For the frequency range 
of interest in the present calculations, 2400-2580 cm -1 , we 
consider the contributions from the v 3 band of C0 2 only and 
exclude any contributions from other bands. From the HIT- 
RAN data base, 18 there are 109 lines listed for this band and 
the maximum value of the rotational quantum numbers 
among them is 108. Although all of them are included in the 


calculations, the diagonization of the anisotropic interaction 
matrix whose size depends on how many states are included 
in the complete set remains tractable. 

The formula used by Boissoles et aL 14 to calculate the 
band-averaged line shape function of v 3 band of C0 2 is simi- 
lar to Eq. (59). However, they introduced an energetically 
corrected function in their formalism. This correction corre- 
sponds to rewriting the spectral density in a more symmetric 
form 19 

/ r o( t * > ) = — Re f e ‘°“ Tr{/*\/p e~‘ L, 4pn}dt. (65) 
TT Jo 

Consequently, the expression of the band-averaged line 
shape function x(oj) is given by 


X(u>) = 4ir 2 n b u'^ l 2 K«lVp/i m l/3)| 2 

a$ m 

X g-^[(C # + C fi )l(G a - Gp) )/2kT 


x r 


R(r) 

dR(r)/dr 


-VUr)lkT\ 


1 



P;lA*,y 


(66) 


Similarly, we follow the same energetically corrected 
procedure to calculate As we mentioned previously, in 
the present calculations we do not treat the velocity average 
exactly. Since the Maxwell ian dist ribution function has a 
sharp maximum at u max = V 2kTlm , we choose 0 a) =u max in 
the trajectory average calculations. With respect to the turn- 
ing points r min , we assume that they are determined by the 
isotropic potential only and that the effect from the aniso- 
tropic potential is averaged out. In the present calculation, 
we simply use the equation, \mv 2 m - V r jso (r inin )=0 to deter- 
mine the turning point r mjn . Of course, more sophisticated 
methods could be used to deal with the velocity averaging 
and the determination of the turning points, however, we do 
not discuss them in the present paper. Thus the formula used 
for the numerical calculations for ?J(io) is 

kT 

&(«) = -8ir 2 n b ij — 2 2 \(o‘\'Ipfi m \/3)\ 2 e flu,2kT 
m „ 

ap m 


T iSO (r) 

1/2 

dR(r) 

1 

[ kT 

Sign 

dr ° a? 


Xr 


1 -2fr( 2+ 4)> v “W 


+ R‘ric; e l 


r = r a^ u '< 



(67) 


As mentioned previously, the value of the parameter 77 intro- 
duced to make the classical trajectory average approximation 
consistent with the closest region approximation is unknown. 
Using a classical model (see Appendix C), we estimate 
77=0.3 which we use in the present calculations. 
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FIG. 1. The band-averaged line shape function x(ai) in units of cm"' atm" 1 
(solid curve) and the correction to the band-averaged line shape function, 
with y=0.3 (dotted curve), as functions of frequency <u (in units of 
cm -1 ) calculated for T- 296 K. 

The calculated band averaged line shape function x(oi) 
and the calculated correction of the band averaged line shape 
function ?<4(oj) are shown in Fig. 1. In this figure we present 
?6{u))I(d instead of since not only has the 

same dimension as x{ui), but also the comparison between 
them makes more sense. The frequency at which '26(<j)/<o 
becomes comparable to x{cd) depends on the short-range an- 
isotropic potential and on the value of the parameter rjr, in the 
present case, this occurs around 45 cm -1 . For smaller fre- 
quencies, ?J(<jj)/<o exceeds ^(oi) implying that the quasistatic 
approximation breaks down and the effects of molecular mo- 
tion must be taken into account. As shown in the Fig. 1, 
?6{u>)lw falls off faster than x(wi) as the frequency increases 
and becomes negligible beyond 80 cm -1 which validates the 
quasistatic approximation for this frequency region. 

The calculated contributions of absorption coefficients 
from the main part F 0 (cu) which corresponds to Boissoles 
et al.’s results 14 are shown in Fig. 2. The total absorption 
coefficients are also shown in Fig. 2 together with the results 
obtained using a Lorentzian line shape and the experimental 
results. 14 In order to investigate how sensitive these results 
are to the choice of the potential parameters, we have re- 
peated the calculations varying the parameters cr, elk, R 2 , 
and A 2 in Eqs. (62) and (63). As expected the results depend 
sensitively on a but less so on elk: small changes in the 
relative weighting of the long- and short-range parts A 2 and 
R 2 can affect the magnitude. We obtained the agreement with 
the experimental data shown in Fig. 3 using the values: 

< 7=3.2 A, e'k = 100 K, R 2 =\.22, and A 2 =0.35. Finally, we 
note that the comparisons discussed above are based on cal- 
culations made with a simplified one-component anisotropic 


FIG. 2. The absorption coefficient o(a») (in units of cm 2 molecule' 1 atm' 1 ) 
vs the frequency at (in units of cm' 1 ) calculated for T =296 K. The curve — 
is the contribution from the band-averaged line shape function and — is the 
total. Also shown are the results calculated using a Lorentzian line shape 
denoted by — and the experimental data denoted by +, respectively. 



FIG. 3. Comparison between the present theoretical results calculated using 
the modified potential parameters and the experimental data. 
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potential. Before accurate quantitative comparisons can be 
made, a more realistic multicomponent potential must be 
used. 

IV. DISCUSSION AND CONCLUSIONS 

As a result of the present work, several conclusions can 
be drawn. First, the effects of molecular motion significantly 
affect the quasistatic line shape for displacements Ao> less 
than 50 cm -1 from the line centers. As noted above for 
smaller displacements, the correction term can exceed the 
quasistatic contributions; however, in this region 
1 /t c <Ao>=1/t where r c and rare time between collision and 
the duration, respectively, both the quasistatic and the impact 
theories break down. As a result, the corrections for molecu- 
lar motion will be more important nearer the band origin or 
in the microwindows between the lines 20 than in the win- 
dows between vibration-rotational bands. Second, the cor- 
rection depend on the interaction potential, and because the 
motion is greater for higher temperatures or for lighter mol- 
ecules, one would expect larger corrections in these cases. 
Conversely, by varying the parameters in the potential in 
order to obtain agreement between theory and experiment, 
one can infer these parameters. However, we note that one of 
the drawbacks of the present quasistatic formalism is that 
only a simplified form of the anisotropic potential having one 
spherical harmonic can be incorporated. Further work to in- 
clude more general potentials is needed before realistic 
potential-energy parameters can be accurately determined. 


Finally, although the present work is a first attempt to con- 
sider the molecular motion a number of simplifying assump- 
tions have been introduced in the calculations; nevertheless 
to take into account of the correction of the quasistatic for- 
malism due to molecular motion is necessary in the near- 
wing region and this becomes a tractable problem especially 
for some simple interacting systems. 

Some refinements of the approximations can be incorpo- 
rated in the formalism. For example, one can improve the 
trajectory average and avoid the introduction of the param- 
eter tj. Because a main uncertainty in the present calculation 
is related to the interaction potential, there are two chal- 
lenges facing us. One is to find a more reliable but still 
tractable potential models to represent the interaction. The 
other is to find a theoretical method to calculate the resolvent 
operator corresponding to more realistic anisotropic poten- 
tials which consists of more than one spherical harmonic 
component. Further work is needed in this direction. 
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APPENDIX A 

The first-order term of the spectral density is 


FJ(o) = j~ f Tr pp ie iU [ 'v\(t')dt',p \dt = - Re [ e im ‘ Tr( pp ie iU \ V[ 

2 irJ-= Jo JJ tr Jo [Jo 


\{f)dt',p 


dt. 


In the step of deriving Eq. (Al), we have used the fact that 


le 


iLl f \\(t')dt',p = —ie iLt f V,(r')*dr',M =ie iU [ v\{-t')*dt' ,p = ie iU f V[(t')dt' ,p 

Jo Jo Jo Jo 


since Vj(r) is an even function as mentioned in the text. Then, 


! ( oj) = Re Tr( p • if V[(t')dt' ,e~ iL ‘pp 1 dt =- Im f Trf p ■ [ ' v\{t')dt' ,e~ iU pp 

If Jo [ [ J 0 jl IT Jo [ [ J 0 


The integration of Eq. (A2) can expressed as follows: 




dt e‘ 


fv'.{t')dt',e~ iU pp [‘v , l (t')dt',e i ^- L) 'pfi\ = r\ f V f l (t , )dt' ,df‘ e i( "- t) ‘' dt' Pt i 
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dt. (A2) 
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J0 [ o 0> — L 

= i(™dt e i0 “ e- iLt V.(t),e' iU —5— p J 
Jo <*>— i, I 

= if”df e i( “- L)t V x (t),-^— dt e‘Wv t (t )— pp, 

Jo a» — L I Jo to — L 


where V'i(f) is a Liouville operator corresponding to V[(r). With Eq. (A3), we rewrite the first-order term as 


Fj(to) = i Re [ e iut Tr fi e~ iL, V x (t) pJdt. 
IT Jo (0 — L 


To further simplify the calculations, we introduce another approximation in Eq. (A4). As discussed in the text, the part of 
^(r), i.e., the term of jt 2 (d 2 r/dt 2 )[dR(r)/dr]G inii0 , which contributes to F,(oj) only, commutes with the time-independent 
interaction V; therefore, e~ lL, V x (t) — e~‘ L °‘V x (t). With this fact and keeping in mind that we expect the effects due to the 
molecular motion to be significant only in near wings, we will approximate the Liouville operator e~ ,Lt by its main part 
e~‘ L °‘ in Eq. (A4). Then we obtain the simple formula 


1 f“ • f - 1 

F x (io)=- Re e ,u “ Tr ft •e"' i o / V' 1 (r) p/x dt 

IT Jo U>~L 


APPENDIX B 

With Eqs. (42) and (49), the trajectory averaged component of the resolvent operator - Z,i)])t, can be 

expressed 

/ i \ /-» v (r) "1 1/2 

Im<iijtl( . T T" 2 r dr{(i x j\){i 2 i\)\ap) 1~~— 

\ W L l / ft.tny. i-Ui a0 J 'min 


I i \ /•» V (r) 1 1/2 

Itn<«ijll( V\ 7 T ) \ i ij\ 1 ) = ~ 2 ' ,rln bVvV2 2 r dr{(i x j\){i 2 i\)\ap) 1--"°— 

x {i"^ r ( 2+r ^ : |( ViS0 ( r ) +/? ( r ) <; ^ ) 


X(^ t |(iij; t )(j2j2))^' £ ^ r /e 4 , (Bl) 

where are the degeneracy factors due to nuclear spin, E t are the energies associated with the state |j), and (2s is the partition 
function for the bath molecules. By carrying out the integration with the help of the S function, we obtain 


MmJiK Vi j 2 (('ijl)(*2‘2)l a 0 t > 

\ 1 / h Trai i-* aB 


X 1- 


'2/2 «P 
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r 
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where r a ^u>') are the positive solutions of the equation u/— /?(r)C a ^=0 which are larger than r mia . In the above expression, 

y. ( r ) r/fR(r) 1 f 1 / \ 

2 <0'iJi)0V2)|a/3 t )' 1_ 1~T Sign r I — rj- I +r drj 


‘2/2 “P 


X(V to (r)+JI(r)C^) 




J. Chem. Phys., Vol. 100, No. 8, 15 April 1994 


5578 Q. Ma and R. H. Tipping: Quasistatic line-shape theory 

With Eqs. (35) and (B2), the correction term of the spectral density, can be rewritten as 

2 If 1 / . 1 \ 1 1 

F 1 (o,)=--Im— da>' Tr ft ■ ■ 3 V; -737- ^7- P (a) M 

7T Z7T{J_ao (it) "a) \ Cl) Z# j ^ ^ t - C li (0 

= ~"2 2 0'l#*l'> ■ 7—— 7 TZjf da)'^ ihrj ,(a)') — — t— — p,.<iW> 

rr ta» o)jj) zirij -00 a>— at a/,'y> 

= "X I 0'l/*l«) ( ~ ~ )3 %,, , y>(at-at,. / ,)p,,(r'|/x|y')- 

ij i* j' ' 

Using the reduced matrix elements. 




1/2 




l(2i+l) 

and reduced functions, £&***! ,(a>), defined by 
1 




(2i+ l)(2i' + 1) 


1/2 _ _ 

2 2 2 2 2 C(Jli;m j mm i )C(j'li',m' j mm' i )^ ijti . j ,(a), 


(B4) 


(B5) 


(B6) 


we complete the summation over the magnetic quantum 
numbers and obtain 


Fl(w) ir 5 ^ (ot-otiy) 3 


</ < 1 


( B? ) 

where the summation indices exclude their magnetic quan- 
tum numbers. For simplicity, since we will use these reduced 
forms instead of the original ones in the results that follow, 
we will omit the superscript redu. 


APPENDIX C 


In general for a specified relative velocity at infinity, 
the trajectory is determined by the impact parameter b and 
the interaction potential V(r), and the relative position of 
molecules along the trajectory s(r) is a function of time 
whose zero has been chosen in the present work as the time 
of closest approach. In order to calculate averages of the 
functions of interest over all trajectories, one has to average 
over b and s(t). However, the functions of interest are usu- 
ally expressed in terms of the intermolecular separation r 
and, therefore, the dynamic variables b and s(t) are not ap- 
propriate and have to be transformed to spherical coordi- 
nates. Assuming that the trajectory is mainly determined by 
the isotropic potential, we have 



00 


2^iso(r) 

mvl 



1/2 


(Cl) 


For a specified r the maximum impact parameter is given by 


b m J.r) = r 


2 V.,Jr) 

1 r~ 


11/2 


mv „ 


(C2) 


With the straight line trajectory approximation, the integral 


of a function F(r) over all trajectories can be expressed as 


f oo C 00 

r=2irj bdbj 
= 4irf 
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F(r)dr 

r»f 2V iso (r) | 1/2 \ 2 j 

= 4 ir 1 j— F(r)r 2 dr, 

J'Jl mv ~ . 


■in 


where r mjn is given by 
2 mv 00~ ^isoC^min) ~ 


b db, 
(C3) 

(C4) 


However, based on the closest region approximation the 
contribution to the integral / comes mainly from the region 
near the turning point of trajectories and, in fact, the function 
F(r) has been expressed in terms of an expansion around 
this point. For consistency, this implies that a cutoff has to be 
enforced on the integration along trajectories. One method to 
introduce the cutoff is to replace 0, the lower limit of the 
integration over b in Eq. (C3), by b min (r) and to assume that 

t> min (r) = b max {r)cos 9 0 , (C5) 

where 0 O is a parameter. Then, one is able to recast Eq. (C3) 
as 


f» fU') 

2 V iSQ (r) b 2 

r F(r)dr 

1 r~ 7 r 

* r min * ^ min ( r ) 

mv a r 


- 1/2 


b db. 



1/2 

F(r)r 2 dr, 


(C6) 


where the parameter 77 finally used in the text is defined by 


77 = sin 0 O = 


1- 



1/2 


(C7) 
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In order to estimate the parameter 77, we consider a 
straight line trajectory with the closest distance denoted by 
rQ ( = 6). in the present case, the function F(r) of interest is 
V[ = ^{d 2 r/dt 2 )(Sf dr) t /, an i S0 ( r )- If we assume that around the 
distance of closest approach the repulsive part of the aniso- 
tropic potential dominates, then F(r)«r -16 and its magni- 
tude drops to 1/2 the maximum value for s=0.3r 0 . With this 
assumption, the interaction takes place along the trajectory 
between -0.3r o <s< + 0.3r o . It corresponds to b tBiB =r 0 
and b m ax = [l + ( s2 / r o)] 1/2r o- Therefore, one is able to esti- 
mate the parameter y=(s/b m3X )—0.3 which we use in the 
calculations. In Eq. (C6) and consequently in Eq. (67) the 
parameter 77 is simply a scale parameter and it is easy to 
calculate results which correspond to different values of rj. 
For example, the correction contribution is comparable to the 
leading contribution at about 27 cm -1 for 77=0.1 while this 
occurs at 45 cm -1 for 77=0.3. Although there is considerable 
uncertainty in the potential for small separations and arbi- 
trariness in the cutoff discussed above and thus in the value 
of 77, we note that the general conclusion that molecular mo- 
tion is important for the near wings of the lines but can be 
ignored in the far-wing region does not depend on the value 
of 77 assumed. 
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